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Abstract 

We report ab initio correlated relativistic calculations of the effective electric field Wd acting 
on the electron in two excited electronic states of PbO, required for extracting the electric dipole 
moment of the electron from an ongoing experiment at Yale, which has the potential of improving 
accuracy for this elusive property by several orders of magnitude. The generalized relativistic 
effective core potential and relativistic coupled cluster methods are used, followed by nonvariational 
one-center restoration of the four-component wavefunction in the heavy atom core. Wd is —3.2 x 
10 24 Hz/(e • cm) for the a(l) state and —9.7 x 10 24 Hz/(e • cm) for the B(l) state. Comparison of 
calculated and experimental values of the hyperfine constant An provides an accuracy check for 
the calculation. 
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Introduction. Following the discovery of the combined CP-parity violation in .Ko-meson 

fl 

decay the search for the electric dipole moment (EDM) of the electron, d e , has become 
one of the most fundamental problems in physics j^J. Considerable experimental effort 
has been invested in measuring atomic EDMs induced by the electron EDM. The best 
available results for the electron EDM were obtained in the atomic Tl experiment , which 
established an upper limit of \d e \ < 1.6 • 10~ 27 e-cm. It is expected that diatomic molecules 
containing a heavy atom can yield more definite results. Modern experiments searching for 
d e in these molecules exploit the fact that the effective electric field seen by an unpaired 
electron, Wd, is greatly enhanced by the relativistic effects relative to the external field 
reachable in a laboratory. The value of Wd is necessary to extract d e from experimental 
measurements. For diatomic molecules with one unpaired electron, such as YbF and BaF, 
semiempirical estimates or ab initio calculations with approximate accounting for correlation 
and relativity provide reasonably reliable Wd values (see Refs. jf], B, ) . These molecules 
are, however, chemical radicals, posing experimental problems. It was pointed out recently 
that the excited a(l) B or B(l) Q states of PbO can be used effectively in the search for 
d e . A novel experiment, using a vapor cell to study excited PbO, has been started at Yale 
University. The unique suitability of PbO for searching the elusive d e is demonstrated by 
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the Yale experiment to the electron EDM, allowing 



four orders of magnitude lower than the current 



the very high statistical sensitivity o 
detection of d P of order 10~ 31 e-cm 

Mt quot ea above. WMe serial calculations R may be va lu ab,e, the authors 
of stressed that "more elaborate calculations were highly desirable". High-order ab 
initio correlated relativistic calculations of the type developed recently 0] are required to 
give accurate values of Wd acting on the unpaired PbO electrons. An accuracy check is 
provided by calculating experimentally known properties which also depend on the electron 
spin density near the heavy nucleus, such as hyperfine constants. 

The terms of interest for PbO in the effective spin-rotational Hamiltonian may be written 

n 

following Ref. [12|. The P,T-odd interaction of d e with Wd is 

H edm = W d d e (J-n), (1) 

where J is the total electron moment and n is the unit vector along the axis from Pb to O. 
The hyperfine interaction of the electrons with the 207 Pb nucleus is 

H Ms = J-A- I, (2) 



where A is the hyperfine tensor, characterized for a linear molecule by the constants A\\ and 
A±, and I is the spin of the 207 Pb nucleus (J = 1/2). 
In practice, the effective operator 



H d = 2d e _ (3) 




is used to express the interaction of d e with the inner molecular electric field E (tr are the 
Pauli matrices), to avoid the large terms which cancel each other because of Schiff's 
theorem. After averaging over the electronic coordinates in the molecular wavefunction, one 
obtains 



h-WTH d (i)\H! a ) , (4) 
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where is the wavefunction for either a(l) or B(l), and Q = (^n\J -nl^n). The hyperfine 
constant A\\ is determined by the expression jlj] 

i 1 

where /ipb is the magnetic moment of 207 Pb, (Xi are the Dirac matrices for the ith electron, 
and T*j is its radius-vector in a coordinate system centered on the Pb atom. 

Both A\\ and Wd depend strongly on the electronic spin density near the heavy nucleus, 
while the molecular bonds are formed in the valence region. As shown previously ( [^l. fllL Il5| 
and references therein) , it is possible to evaluate the electronic wavefunction near the heavy 
nucleus in two steps. Using this strategy here, a high-accuracy relativistic coupled cluster 
(RCC) calculation of the molecular electronic structure with the generalized relativistic 
effective core potential (GRECP) is carried out, providing proper electronic density in the 
valence and outer core regions. This is followed by restoration of the proper shape of the 
four-component molecular spinors in the inner core region of the heavy atom. 

Methods and calculations. A 22-electron GRECP for Pb [17] is used in the first stage of 
the two-step calculations of PbO: the inner shells of the Pb atom (Is to 4/) are absorbed 
into the GRECP, and the 5s5p5d6s6p electrons and all the oxygen electrons are treated 
explicitly. Two series of calculations are carried out, denoted as (a) and (b): calculation (a) 
correlates 10 electrons, freezing the 5s5p5d shells of Pb and the Is shell of O; (b) correlates all 
30 electrons treated explicitly. States with the leading configurations afa^Tif, afa^ufnl, an d 
°"i cr 2 7r i 7r 2 are calculated. Here <7i )2 and 7r 12 are molecular valence orbitals, with the subscript 



enumerating them in order of increasing energy. For each series of calculations, correlation 
spin-orbital basis sets are optimized in atomic two-component GRECP/RCC calculations 
of Pb. The four 6s and 6p electrons are correlated in the basis set optimization stage of 
calculation (a), and 22 electrons (5s to 6p) are correlated in the optimization of the basis 
set used in series (b). Correlation is taken into account at this stage by the RCC method 
with single and double excitations (RCC-SD) jlgj ; the average energy of the five lowest 
states of Pb is minimized. The detailed description of the basis set generation procedure 



may be found in Refs. 
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20(. A [4s3p2<i] basis, obtained by omitting the / function from 



Dunnin g's correlation-consistent (10sbp2dlf)/[4s3p2dlf] basis listed in the MOLCAS 4.1 



library 



2l|, is used for oxygen. We found that the / orbital has little effect on the core 



properties calculated here. Previous calculations show that these basis sets are adequate for 
our purpose. 

PbO calculations start with a one-component self consistent field (SCF) computation 
of the molecular ground state, using the spin-averaged GRECP (AGREP). The Pb spinors 
5s5p5d are frozen in the (a) series, using the level-shift technique |22|. An AGREP/RASSCF 



(restricted active space SCF) calculation 0, Q of the lowest S + state of PbO is then 
performed. In the RASSCF method, orbitals are divided into three active subspaces: RAS1, 
with a restricted number of holes allowed; RAS2, where all possible occupations are included; 
and RAS3, with an upper limit on the number of electrons. 

Different distributions of electrons in these active subspaces are used (details may be 
found in j^) to estimate the different correlation contributions to the RASSCF values of 
An and Wa- Two-component RCC-SD molecular calculations are then performed. The 
AGREP/RASSCF calculations include only the most important correlation and scalar- 
relativistic effects, while the GRECP/RCC-SD calculations also account for spin-orbit in- 
teraction. The Fock-space RCC calculations start from the ground state of PbO and use 
the scheme 

PbO + «- PbO -> PbO" 

\ / (6) 

PbO* 



Details on the model spaces used may be found in |24|. 

Only valence and outer core electrons have been treated so far. Since we are interested in 
properties near the Pb nucleus, the shape of the four-component molecular spinors has to be 
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restored in the inner core region. Al 
one-center restoration scheme (see 
two steps: 



molecular spinors are restored using the nonvariational 



22|, |25[ and references therein). This is done in 



First, equivalent numerical one-center basis sets of four-component spinors and two- 
component pseudospinors are generated by the finite-difference all-electron Dirac-Hartree- 
Fock (DHF) and GRECP/SCF calculations, respectively, of the same valence configurations 
of Pb and its ions. In the DHF calculations the inner core spinors (Is to 4f) are frozen after 
the calculation of Pb 2+ , and the nucleus is modeled by a uniform charge distribution within 
a sphere of radius r nuc i = 7.12fm = 1.35 x 10 _4 a.u. The root mean square radius of the 
nucleus is 5.52 fm, in accord with the parametrization of Johnson and Soff [2g|, and agrees 
with the 207 Pb nucleus experimental value of 5.497 fm 2?|. Taking the experimental value 
for the root mean square radius and a Fermi distribution for the nuclear charge changes 



An and W d by 0. 



GRECP/HFJ 17 
veloped in Refs 



or less. The all-electron four-component HFD [28] and two-component 



291 codes are employed for the basis generation, using the procedure de- 
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The basis sets generated are [9sl4p7ci] for series (a) and [6s7p5d] 
for series (b), with the latter carefully optimized. These sets are orthogonal to the inner 
core (see above). They describe mainly the core region, and are generated independently of 
the basis set for the molecular GRECP calculations discussed earlier. 

In the second step, the basis of one-center two-component atomic pseudospinors is used 
to expand the molecular pseudospinorbitals; these two-component pseudospinors are then 
replaced by the equivalent four-component spinors, retaining the expansion coefficients. A 
very good description of the wave function in the core region is obtained. 
The RCC-SD calculation of Wd and An employs the finite field method 



3fl|3l|. The 



operator corresponding to the desired property [Eq. (0) or (J2J)] is multiplied by a parameter 
A and added to the Hamiltonian. The first derivative of the calculated energy with respect 
to A gives the evaluated property. This is strictly correct only at the limit of vanishing A, 
but it is usually possible to find a range of A where the energy is linear in A and the energy 
changes are large enough to attain the required precision. The quadratic dependence of the 
energy on A is eliminated here by averaging the components of a given term, a(l) or B(l), 
with opposite signs of A. 

Results and discussion. Calculated results for the (a) and (b) series are presented in 
Tabled The internuclear distance is 2.0 A. The RASSCF calculations use the 22-electron 
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TABLE I: Calculated parameters A\\ (in MHz) and W d (in 10 24 Hz/(e • cm)) for the o(l) and B(l) 
states of 207 PbO. The experimental value of An in a(l) is —4113MHz. The preliminary value of 
An in B{1) is 5000 ± 200MHz 1 

State a(l) ofofTif^ 3 Xi B(l) o\a\^\ 3 U X 

Parameters A\\ Wd An Wd 



Expansion s s,p s,p,d s,p s,p,d s s,p s,p,d s,p s,p,d 

lOe-RASSCF -759 -1705 -1699 0.96 0.91 1900 0.0 0.0 

lOe-RCC-SD -2635 -2.93 3878 -11.1 

30e-RCC-SD -359 -3062 -3012 -3.08 -3.18 195 4510 4568 -10.4 -9.7 
1 private communication, D.DeMille, 2003. 

GRECP for Pb. Twenty of the 30 electrons treated were in the inactive space, and only 10 
were correlated. Using the C2 V classification scheme, 2 Ai orbitals are in RAS1, 6 orbitals 
(2 Ai, 2 Bi, and 2 B 2 ) in RAS2, and 41 (16 A 1; 5 A 2 , 10 B 1; and 10 B 2 ) in RAS3. No 
more than two holes in RAS1 and two particles in RAS3 are allowed. The basis sets on Pb 
are (14sl8pl6d8/)/[4s7p5d3/] for the RASSCF and 30-electron RCC-SD calculations and 
(15sl^l2d9/)/[5s7p4d2/] for 10-electron RCC-SD. A (10s5p2d) /[As3p2d] basis is put on O 
in all calculations. 

We discuss mainly the results for the a(l) state (leading configuration o-fo-^nfiil), for 



which the reliable experimental value of An is available (—4113 MHz) 32j and a semiempirical 
estimate of \Wd\ > 12 x 10 24 Hz/(e-cm) was made recently Q]. There are several points to 
note: (1) Inclusion of the spin-orbit interaction changes A\\ and Wd dramatically, as may be 
seen from the difference between the 10-electron RASSCF and RCC-SD results. (2) The ab 
initio value of Wd is four times smaller than the semiempirical estimate ^(|. (3) Accounting 
for outer core-valence correlation by 30-electron RCC-SD changes Wd by 5% and A\\ by 
15%, yet the error in the calculated An is 25%; calculations on BaF |j| and YbF [8] gave 
10% accuracy. (4) An is mainly determined by the p wave, whereas Wd mostly comes from 
s-p mixing. 

The need for including correlation in the PbO molecule for the properties discussed here 
can be seen already in the semiempirical model Q]. The leading contribution to the highest 
occupied 02 orbital in this model comes from the Pb 6s atomic orbital, with a weight of ~0.5 
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(the corresponding coefficient in the molecular orbital expanded as a linear combination of 
atomic orbitals is ~0.7). This contradicts the qualitative analysis of the chemical bond 
formation, which predicts that the Oi orbital is mainly formed from the oxygen 2p a and lead 
6p a orbitals. The RASSCF calculations of the lowest 3 S + state confirm this point, with the 
weight of the Pb 6s orbital varying between 0.04 for 10 active electrons and 0.1 for 30 active 
electrons. The weight of the oxygen 2p a is ~0.5 and that of the lead 6p a is ~0.1, whereas 
o~i consists mainly of the lead 6s orbital, with negligible contribution from lead 6p a . Note 
that the oxygen 2p a and lead 6p a orbitals are not ortho gon al to each other; after one-center 
reexpansion of the oxygen basis functions on lead (see and Eq. (6) in Hj), the weight 
of the 6p a orbital goes up to 0.3. We expect that such strong admixture of the s-wave to 
the o"2 orbital would not appear in the semiempirical model if configurations describing the 
correlation of the o-i electrons were included in the model space. It is important to add that 
the lowest virtual 03 orbital gets the main contribution from the lead 6p a (with a weight of 
about 0.5), and the configurations containing this orbital are first admixed into the leading 
configuration of the a(l) state due to the spin-orbit interaction on Pb. 

If the spin-orbit interaction is neglected, the s-wave contribution to An and the s,p- 
wave contributions to W4 are due primarily to correlation of the a electrons. The RASSCF 
calculation indicates (see Table HJ) that such contributions increase A\\ but decrease Wd, 
resulting in a sign change for Wd, in agreement with the final RCC-SD result (details may 
be found in Ref. j^j])- Besides, as correlation is expected to have a strong influence on 
the values of A\\ and Wd, introducing the SO interaction with the 3 II and 1 II states by just 
mixing the corresponding a and 7r orbitals may not be satisfactory. All these conclusions 
could be reached only after extensive molecular calculations, and the estimates made in 
Ref. were important at the first stage of the experimental effort. 

As may be expected, the accuracy of the calculated A\\ and Wd values is lower for such a 
complicated system as the excited states of the PbO molecule than for the ground states of 
BaF and YbF. The valence electron in the latter molecules is in a a orbital, with much higher 
density near the heavy nucleus than the valence ir electrons in PbO. Thus, the s,p, d-waves 
on the Pb nucleus are affected more strongly by correlation, and higher-order inclusion of 
correlation (triple and quadruple amplitudes in the RCC method) as well as larger basis 
sets may be necessary. As pointed out above, Wd is more stable than An with respect to 
changing the number of correlated electrons, and we expect the accuracy of the calculated 



7 



Wd to be better than for A\\. Our estimated error bounds put the real Wd between 75% 
and 150% of the calculated value, which is quite satisfactory for the first stage of the EDM 
experiment on PbO. It should be noted that the estimate of statistical sensitivity to the 
electron EDM made in is based on a Wd value close to that obtained here. 

A detailed analysis of correlation and spin-orbital effects on A\\ and Wd in PbO will be 
published elsewhere |33[. Unfortunately, the experimentally available Aii of the a(l) and B{\) 
states provides a check on the p wave only. It would be desirable to measure Ah in some state 
with an excited <7i-electron, with the main contribution coming from the s-wave. Another 
accuracy check, using A\\A±, is not applicable here, because of experimental difficulties in 
measuring the very small A± for diatomic molecules with total electronic momentum J > 1. 
Our estimate of the accuracy of the calculated Wd is therefore not as straightforward as for 
YbF and BaF Q, j). 

Finally, we would like to note that we identified the lowest 3 IIi state as B(l) according 
to the AS classification given in 3^. Conclusive identification requires more extensive ab 
initio correlation calculations. 
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